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ABSTRACT 

We present comparison of X-ray proxies for the total cluster mass, including the spectral temperature (Tx), 
gas mass measured within rsoo iM„), and the new proxy, Yx, which is a simple product of Tx and M„ and is 
related to the total thermal energy of the ICM. We use mock Chandra images constructed for a sample of clus- 
ters simulated with the eulerian A^-body-ngasdynamics adaptive mesh refinement ART code in the concordance 
ACDM cosmology. The simulations achieve high spatial and mass resolution and include radiative cooling, 
star formation, and other processes accompanying galaxy formation. Our analysis shows that simulated clusters 
exhibit a high degree of regularity and tight correlations between the considered observables and total mass. 
The normalizations of the M - Tx, M„ - Tx, and M - Yx relations agree to better than 10 - 15% with the 
current observational measurements of these relations. Our results show that Yx is the best mass proxy with a 
remarkably low scatter of only » 5 - 7% in M500 for a fixed Yx, at both low and high redshifts and regardless 
of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx - M500 relation is 
close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the 
cluster mass function from X-ray selected samples. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: clusters - clusters: formation - methods: 
numerical 



1. INTRODUCTION 

The evolution of the cluster abundance is one of the most 
sensitive probes of cosmology, which can constrain the power 
spectrum normalization, matter content, and the equation of 
state of the dark energy. The potential and importance of 
these constraints have motivated efforts to construct several 
large surveys of high-redshift clusters during the next several 
years. However, in order to realize the full statistical power 
of the upcoming cluster surveys, it is paramount that the rela- 
tion between cluster mass and observables and any potential 
biases are well known. 

Several cluster observables based on the galaxy velocities, 
optical light. X-ray observables such as luminosity, temper- 
ature, mass of the intracluster medium (ICM), and Sunyaev- 
Zel'dovich (SZ) flux have been proposed and used in the lit- 
erature as proxies of the total cluster mass (see Voit 2005, for 
a recent comprehensive review). In this study we focus on 
the mass indicators derived from cluster X-ray observables, 
which provide a handle on the properties of the hot ICM com- 
ponent. X-ray luminosity, Lx, computed using the flux inte- 
grated within a certain radius or a range of radii, is expected to 
correlate with cluster mass (e.g.. Kaiser 1986 ) and is the most 
straightforward mass indicator to measure observationally. Lx 
has been used for cosmoiogicai fits to the cluster samples 
from ROSAT All -Sky Survey (Reiprich & Bohringer 200^ 
[ Allen et all 120031) and Deep Cluster Survey (Borgani et al. 
I2001h . However, Lx is also the least accurate (internally) 
of all proposed X-ray proxies for Mtot. Lx is dominated by 
the cluster cores and thus is particularly susceptible to non- 
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gravitational processes i n flie IC M. Given the large sca tter in 
the Lx-Tx relation (e.g. jDavid et al. 1993; Markevitch ri998t 
llkebe e t al. 2002), the Lx- M relation for real clusters proba- 
bly also has significant scatter (.Stanek et al.^2006) . The slope 
of the Lx - M relation dev iates from the self-similar predic- 
tion (e.g., 'Allen et al.'"2003"). In addition. X-ray luminosity is 
notoriously difficult to reliably model in cosmoiogicai simu- 
lations (e.g., Anninos & Norman 1996; Lewis et al. 2000), a 
significant disadvantage given that simulations are often used 
to get a handle on the expected evolution of the mass vs. proxy 
relations. These problems could potentially be alleviated with 
sufficient angular resolution by excising t he flux f rom cluster 
cores, responsible for most of the scatter (Markevitch 1998). 

The most common choice of mass proxy used to mea- 
sure the cluster number density and constrain cosmoiogi- 
cai parameters is the X- ray temperature of the intraclus- 
ter plasma (e.g. .[Henry & Arnaud 1991; Oukbir & Blanchar3 
1992; Markevitch 1998; Henry 2000; Seljak 2002; Ikeb e et all 
.2002; Pie rpaoli et al. 2003). Until recently, there was a 
large apparent systematic uncertainty in the normalization 
of the M - Tx relation, as evidenced, for example, by a 
~ 30 - 50% discrepancy between obser vational measure- 
ments and cosmoiogicai simulations (e.g.. Finogueno v et alJ 
i2001; Pierpaoli et al. 200 3). Over the last several years the 
M - Tx normalization was revised both in simulations and 
observations due to (1) inclusion of more realistic physics 
in cosmoiog icai simulations (e.g., radiative co oling and star 
formation, Dave et al. 2002, Muanwong et al. 2002), (2) im- 
proved analyses of observed clusters using more realistic 
gas d ensity profiles (e.g., Borgani et al. 2004; Vikhlinin et al] 
mMi . (3)' more reliable measurements of the cluster temper- 
ature profile s (Markevitch et al. 1998; Nevalainen et al. 200(J 
'Arnaud et al. 2005; VikhHnin et al. 2006), and (4) better un- 
derstanding of the meaning of the mean spectral X-ray tem- 
perature, Tx, and the use of uniform de finition of Tx in ob- 
servations and in simulations analyses ( Mazzotta et alJl2()04l 
Rasia et al. 2005; Vikhlinin 2006). The current agreement be- 
tween models and observations is « 10% (see below). 
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The scatter in the M - Tx relation is significantly reduced 
compared to that in the Lx- M relation (the upper limit from 
observations is w 15% in M for fixed T for relaxed clusters; 
[Vilchlinin et al. 2006). In general, existence of a tight relation 
such as M - Tx indicates that clusters are remarkably regular 
population of objects with their global properties tightly re- 
lated to total mass (e.g., Mohr et al. 1999), and scatter caused 
by secondary effects such as substructure in the ICM, non- 
gravitational processes, and mergers (O'Hara et al. 2006). 

More recently, gas m ass was u sed as a proxy for Mtot 
dVikhlin in et al. 2003; Voev odkin & V ikhlinin 2004). The 
practical advantage of gas mass over temperature is that it can 
be measured robustly from the X-ray imaging alone. Also, the 
Mg - M relation in principle permits external calibration from 
the CMB measurements of the global baryon-to-dark matter 
ratio. Finally, Mg can be expected to be less sensitive to merg- 
ers which should translate into smaller scatter in the Mg - M 
relation. The caveat is that trend of gas mass with cluster mass 
and evolution with redshift are not yet fully understood. 

The use of clusters as efficient probes for precision cos- 
mology puts stringent requirements on observable cluster 
mass proxies: 1) tight, low-scatter correlation between the 
proxy and mass, 2) with the scatter insensitive to merg- 
ers, as the frequency of mergers is expected to increas e 
shai-ply with redshift (iGottlober. Klvpin. & Kravtsovll2001h . 
3) simple power-law relation and evolution which can be de- 
scribed by a small number of parameters and be as close 
as possible to the prediction of the self-similar model. The 
point (3) is crucial to ensure that the self-calibration strate- 

f'es for analyses of large cluster surveys ("Levine et al. 2002^ 
u 2003; Maiumdar & Mohr 2003, 2004; Lima & Hu 2004, 
)05; Wan2 et al. 2004) are successful. It is particularly im- 
portant that the sc atter in the observa ble-mass relation is low 
and well-behaved jLima & HJ2005h . 

In general, a mass proxy does not have to be a single clus- 
ter property, such as Lx, Tx or Mg. Any physically-motivated 
combination of these variables that is expected to be tightly 
related to cluster mass can be used to construct a valid mass 
indicator. A hint for a better X-ray mass proxy is provided 
by recent studies based on cosmological simulations of clus- 
ter formation ( Motl et al. 2005; Nagai 2006), which show that 
integrated SZ flux, Ysz, proportional to the product of gas 
mass and temperature, is a good, robust mass indicator with 
low scatter in the Ysz - M relation, regardless of the dynam- 
ical state of the cluster. In addition, the Ysz - M relation 
exhibits a simple, nearly se lf-similar evolution with redshift 
ida Sil va et al. 2004: Nagai 2006). The physical reason for 
the robustness of the SZ flux is straightforward: Ysz is directly 
related to the total thermal energy of the ICM and thus to the 
depth of the cluster potential well (see eq.0]in §|2lbelow). 

In this study, we show that a similar robust, low-scatter 
mass indicator can be constructed using X-ray observables. 
The indicator, which is simply the product of the total ICM 
mass and X-ray spectroscopic temperature, Yx = MgTx, cor- 
relates strongly with cluster mass with only w 5 -8% intrinsic 
scatter. The scatter is robust to mergers, in the sense that even 
for disturbed unrelaxed systems it gives unbiased estimates of 
mass with the statistical uncertainty similar to that for relaxed 
systems. Thus, the scatter of the Yx - M relation at higher 
redshift is similar to the scatter at z = 0. In addition, we show 
that evolution of the slope and normalization of the Yx-M re- 
lation is nearly self-similar. These properties make Yx partic- 
ularly useful for measurements of cluster mass function using 
X-ray surveys. 



2. MASS PROXIES 

Physical properties of virialized systems, such as clusters, 
are expected to correlate with their total mass. For example, 
in the self-similar model ( Kaiser.,. 1986, 1991) the cluster gas 
mass is expected to be simply proportional to the total mass: 

Ma, = CM.Mgas, (1) 
where masses are determined within a radius enclosing a cer- 
tain overdensity Ac with respect to the critical density of the 
universe at the epoch of observation, Pmt{z), and Cm, is a 
constant independent of cluster mass and redshift. The self- 
similar relation between cluster mass and temperature is: 

E(z)Ma^.^CtT''\ (2) 

Here the function E{z) = H{z)/Hq for a flat cosmology with 
the cosmologic al constant a ssumed throughout this study is 
given by (e.g.. lPeeblesll993h : 

E(z) = ^JnM(l+zy+^A, (3) 
where and Qa are the present-day density parameters for 
matter and cosmological constant. 

The SZ flux integrated within a certain radius, Ysz, is pro- 
portional to the total thermal energy of the ICM gas and thus 
to the overall cluster potential, which makes it relatively in- 
sensitive to the details of the ICM physics and merging: 

Ysz^l'^] r nJ,dV^MgT„„ (4) 

\ nieC^ j Jv 

where ks, crj, nie, and c have their usual meaning, are Tg 
the electron number density and temperature of the gas, and 
T„, is the gas mass-weighted mean temperature of the ICM. 
Combinations of equations J1I2I4> gives self-similar predic- 
tion for the Fsz - M relation, 

£(z)2/5Ma, - Cy,jI^' (5) 
Cosmological simulations show that Ysz is a good, low- 
scatter cluster mass proxy and that Ysz - M relation 
form and evolution are close to the self-similar prediction 
(da Silva etal. 2004; Motl et al. 2005; Hallmanetal. 2006| 
Nagai 2006). Given the good qualities of Ysz as a mass proxy, 
it is interesting whether a similar indicator can be constructed 
from the X-ray observables, which could be used in studies of 
the X-ray cluster abundances. The simplest X-ray analog of 
Ysz is 

Yx = Mg Tx, (6) 

where Mg is the gas mass derived from the X-ray imaging data 
(it is measured within a radius enclosing overdensity A^) , and 
Tx is the mean X-ray spectral temperature (Mazzott a et alJ 
I2004t lVikhlininli2006.) . As we describe below in § it is 
advantageous to measure Tx excluding the central cluster re- 
gion, which can be achieved with moderate angular resolu- 
tion X-ray telescopes (< 15" FWHM). To excise the central 
regions is desirable because the observed cluster temperature 
profiles show a greater degree of similarity outside the core 
(Vikhlinin et al. 2006), and also because this makes the spec- 
tral temperature closer to the gas mass averaged r,„ which 
ideally should be used in equation|6l 

3. MOCK CHANDRA IMAGES AND ANALYSES 
OF SIMULATED CLUSTERS 

A detailed account of the simulations, mock image gen- 
eration, and analysis will be presented elsewhere (Nagai, 
Vikhlinin & Kravtsov 2006, in preparation). Here, we give 
a brief overview of the procedure and define how the observ- 
ables used in this study are derived. 
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TABLE 1 

Simulated cluster sample at z = 



Cluster 
ID 


500 

(10"/i" 


Mil 
500 

' Mo) 


(Tx) 
(keV) 


CLIOI 


90 8 


8 17 


0. / 


CL102 


54 5 


4.82 


5 8 


CL103 


57 1 


4 91 


4 8 


CL104 


53 9 


5 15 


7 7 


CL105 


48.6 


4.71 


6.2 


CL106 


34.7 


3.17 


4.3 


CL107 


25.7 


2.17 


3.9 


CL3 


20.9 


1.91 


3.6 


CL5 


13.1 


1.06 


2.4 


CL6 


16.8 


1.38 


3.4 


CL7 


14.1 


1.21 


2.9 


CL9 


8.23 


0.73 


1.6 


CLIO 


6.72 


0.43 


1.9 


CLll 


8.99 


0.78 


2.0 


CL14 


7.69 


0.62 


1.8 


CL24 


3.47 


0.26 


LO 



3.1. Simulated Cluster Sample 

In this study, we use high-resolution cosmological simu- 
lations of sixteen cluster-sized systems in the "concordance" 
flat ACDM model: = 1 - Qa = 0.3, Qb = 0.04286, 
h - 0.7 and erg - 0.9, where the Hubble constant is defined 
as lOO/i km s"' Mpc"', and erg is the power spectrum normal- 
ization on 8/j ' Mpc scale. The simulations were done with 
the A daptive Refinernent Tree (ART) A^-bo dy-Hgasdynamics 
code ('Kravtsov< ll999l iKravtsov et alJl200a) . a Eulerian code 
that uses adaptive refinement in space and time, and (non- 
adaptive) refinement in mass to reach the high dynamic range 
required to resolve cores of halos formed in self-consistent 
cosmological simulations. 

To set up initial conditions we first ran low resolution sim- 
ulations of two 80 /i"' Mpc boxes and seven 120/;"' Mpc 
boxes, from which we selected sixteen clusters with the virial 
masses ranging from « 7 x lO'"* to 2 x lO'*" Mq for 
re-simulation at higher resolution. High-resolution simula- 
tions were run using 128^ uniform grid and 8 levels of mesh 
refinement in the computational boxes of 120 /i"' Mpc for 
CLlOl-107 and 80/1-' Mpc for CL3-24. These simulations 
achieve the dynamic range of 128 x 2'^ = 32768 and peak 
formal resolution of ^ 3.66 h^^ kpc and 2.44 h'^ kpc, corre- 
sponding to the actual resolution of 7/; ' kpc 5/i ' kpc 
for 120 and 80/? ' Mpc boxes, respectively. Only the re- 
gion of ~ (3 - 10) Mpc around the cluster was adap- 
tively refined, the rest of the volume was followed on the 
uniform 128^' grid. The mass resolution corresponds to the 
effective 512^ particles in the entire box, or the Nyquist wave- 
length /iNy = 0.469 /z"' and 0.312/? ' comoving megaparsec 
for CLlOl-107 and CL3-24, respectively, or 0.018/1' and 
0.006 /z ' Mpc in the physical units at the initial redshift of 
the simulations. The dark matter particle mass in the region 
around the cluster was 9.1 x 10**/; ' Mq for CLlOl-107 and 
2.7 X 10** /z"' Mq for CL3-24, while other regions were simu- 
lated with lower mass resolution. 

The cluster simulations used in this analysis include dissi- 
pationless dynamics of dark matter, gasdynamics, star forma- 
tion, metal enrichment and thermal feedback due to the su- 
pernovae type II and type la, self-consistent advection of met- 
als, metallicity dependent radiative cooling and UV heating 
due t o cosmological ionizing background iHaardt & Madau. 
Il996h . The cooling and heating rates take into account Comp- 



ton heating and cooling of plasma, UV heating, atomic and 
molecular cooling and are tabulated for the temperature range 
10^ < T < 10^ K and a grid of metallicities, and UV inten- 
sities using the cloudy code (ver 96b4, Ferland et al. 1998). 
The CLOUDY cooling and heating rates take into account metal- 
licity of the gas, which is calculated self-consistently in the 
simulation, so that the local cooling rates depend on the local 
metallicity of the gas. 

Star formation in these simulations was i mplemented usin g 
the observationally-motivated recipe (e.g., Kennicutt 1998"): 
P* - P^-ai ^* with f» = 4 X 10^ yrs. Stars are allowed to form 
in regions with temperature T < 2 x lO'^K and gas density 
n > 0.1 cm"^. No other criteria (like the collapse condition 
V ■ V < 0) are used. Comparison of the runs with different 
choices of the threshold for star formation, n - 10, 1, 0.1, 
and 0.01 cm""*, shows that the threshold affects gas fractions 
at small radii, r/ryir < 0.1, but its effect is negligible at the 
radii we consider in this study. Further details about this sim- 
ulated cluster sample will be presented in Nagai, Vikhlinin & 
Kravtsov (2006, in preparation). 

The properties of the simulated cluster sample are summa- 
rized in Table[n We list the total and hot ICM masses defined 
within a radius rsm enclosing the overdensity of 500 with re- 
spect to the critical density at z = 0. We also list the mean 
spectral temperature measured using mock Chandra spectra 
from the radial region of (0.15 - 1) rsoo, as described in § 13.31 
Although the simulated cluster sample is small, the objects 
cover more than an order of magnitude in mass and are sim- 
ulated with very high resolution which allows us to take into 
account effects of galaxy formation on the ICM. 

In what follows, we use cluster total mass and observables 
measured within rgoo ~ 0.5-0.6 ryir, where is a traditional 
definition of cluster virial radius using the virial overdensity 
Avii- (~ 337 at z = for the cosmology adopted in our sim- 
ulations) with respect to the mean density at z = 0. This 
choice of the outer radius is mainly motivated by the fact that 
cluste rs are more relaxed within rgpn comp ared to the outer re- 
gions dEvrard et alJll996tlLau et alJ2006 ). Also, the analysis 
of real clusters is often limited to a similar radius because of 
the limited field of view. 

3.2. Mock Chandra Images 

We first create the X-ray flux maps of the simulated clus- 
ters viewed along three orthogonal projections. The flux map 
is computed by projecting the X-ray emission arising from the 
computational cells within 3 /?vir of the cluster along a given 
line-of-sight, taking into account the actual gas density, tem- 
perature, and metallicity of each cell in a simulation output. 
We compute the X-ray plasma emission, A£(r/,Z;, z), using 
the MEKAL code with the relative elemental abundance from 
Anders & Grevesse ( 1989), and assuming interstellar absorp- 
tion with a hydrogen column density ofnn = 2.0 x 10^" cm"^. 
This provides expected emission spectra in the 0.1-10 keV 
energy range on a spatial grid 1024 x 1024 pixels with the 
size 4.88 h-^ kpc for CLlOl-107 and 2.44 h'^ kpc for CL3- 
24. The entire map is therefore 5.0 /i"' Mpc and 2.5 /i"' Mpc, 
respectively. In generating the spectral maps we assume the 
redshift of cluster observation is Zobs - 0.06 for the simulation 
output at z = and Zobs = 0.6 for the z = 0.6 sample. 

Next, we simulate the photon map by convolving the spec- 
trum from each image pixel with the on-axis response of the 
Chandra ACIS-I camera and drawing the number of photons 
in each spectral channel from the Poisson distribution. We 
simulate two sets of the mock Chandra photon maps. In the 



4 



KRAVTSOV, VIKHLININ & NAGAI 



first set we assume a 100 ksec Chandra exposure, which is 
fairly typical for the real deep observations. From this set, 
we generate cluster images in the 0.7-2 keV band. We then 
add the uniform Pois son background with t he intensity typ- 
ical for ACIS-I data (Markevitch et al. 2003). These images 
are used to identify and mask out from the further analysis all 
detectable small-scale clumps, as routinely done by observers. 
These clumps contain a small fraction of gas mass and do not 
bias gas mass estimates. They can, however, bias tempera- 
ture measurement significantly. Our clump detection is fully 
automatic and based on the wavelet decomposition algorithm 
described in Vikhlin in et al.l lfT998). 

All further analysis is performed using the second set of 
photon maps generated for very long exposures, 10'' sec for 
the z - Q sample and 10^ sec exposure for the z - 0.6 sam- 
ple. Here, the exposure time is chosen to ensure ^ 10'' - 10^ 
photons outside the cluster core region for all simulated clus- 
ters. This second set of data is used in our analysis to derive 
gas mass and mean spectroscopic temperature of the ICM. 
The exposures are artificially large by design as we are in- 
terested in the intrinsic scatter of the X-ray observable mass 
relation, not the statistical errors due photon noise in a partic- 
ular choice of short exposure. Also, we ignore further com- 
plications present in reduction of real Chandra data, including 
background subtraction and spatial variations of the effective 
area (i.e., we assume that accurate corrections for these effects 
can be applied to the real data and any associated uncertainties 
are included in the reported measurement errors). 

3.3. Analysis of Mock Chandra Data 

We analyze the mock data for the simulated clusters em- 
ploying the techniques and tools used in analyses of the 
real Chandra cluster observations, as described in detail by 
[Vikhlinin et al. (2006). After masking out detectable sub- 
structures (see above), we fit the X-ray spectra in concentric 
annuli to measure the projected temperature and metallicity 
profiles. Next, we measure the X-ray surface brightness pro- 
file in the 0.7-2 keV band (used in the real data analysis be- 
cause it provides the optimal signal-to-noise ratio). Using the 
effective area as a function of energy and observed projected 
temperature and metallicity at each radius, we convert the ob- 
served brightness from units of count rate to the projected 
emission measure, EM - J ngnp dl. The derived projected 
emission measure profile is fit to a three-dimensional model 
of Pgas('") that has great functional freedom and can indepen- 
dently describe the gas density slopes at both r > rgoo, and at 
smaller radii and in the very inner region. The best fit model 
directly provides the gas mass profile, Mgas(r). 

We also measure the average X-ray spectral temperature, 
Tx from a single-temperature fit to the spectrum integrated 
within r5oo, excluding the central region strongly affected by 
radiative cooling. The inner cut is set at a fixed fraction rsoo, 
Tin = 0.15 r5oo. Note that we choose to cut out central region 
defined using a fixed f raction of r^on rather th an a fixed metric 
radius of 70 kpc as in IVikhlinin etar('2()()6^. This new defi- 
nition for Tx results in only a small correcti on (- 3% on aver- 
age) of the Tx values reported in Vikhl inin et al.l f2006. who 
used notation of Tspec instead of Tx)- We choose the specific 
value of Tin = 0.15 r5oo because beyond this radius the ob- 
served temperature profiles of clusters are approximately self- 
similar, while at smaller radii the profiles show a large scatter. 
The choice of ri„ thus should maximize the self-similarity of 
the relation between mass and Tx- 

In our analysis below we distinguish the unrelaxed clusters 



from the relaxed systems to test for sensitivity of the mass 
proxies to mergers and substructure. Specifically, as is usually 
done to classify observed clusters, we visually examine mock 
100 ksec Chandra images for x, y, and z projections of each 
cluster and classify as "relaxed" clusters which have regular 
X-ray morphology and no secondary maxima and minimal de- 
viations from the elliptical symmetry. "Unrelaxed" clusters 
are those with secondary maxima, filamentary X-ray struc- 
tures, or significant isophotal centroid shifts. 

4. COMPARISON OF MASS INDICATORS 

In this section we compare three indicators of cluster mass: 
X-ray spectral temperature Tx, gas mass Mg, integrated SZ 
flux Ysz, and its X-ray analog Fx = Mg Tx- 

Figures[lH3lshow the M500 - Mg, M500 - Tx, and M500 - Fx 
relations, respectively. We present correlations at redshifts 
z = 0.0 and z - 0.6, classifying clusters into relaxed and un- 
relaxed based on their X-ray image morphology, as described 
in § 13.31 We performed power law fits to these relations and 
present best fit values of parameters, as well as the amount 
of scatter around the best fit relation for different subsets of 
clusters in Table |2] The upper half of the table shows fits for 
best fit power law normalization and scatter for the values of 
slope fixed to the value predicted by self-similar model. The 
lower half shows best fit values for the case when both nor- 
malization and slopes were allowed to vary. 

For comparison, we also include fits to the Fsz - M rela- 
tion, where Fsz is th e three-dime nsional integrated SZ flux 
measured within rsoo (*N agaill2'006h . The power law fits were 
performed for each projection (i.e., x, y, and z) separately. We 
estimate the uncertainties in the best fit parameters by gener- 
ating 10000 bootstrap samples CPress et al. 1992) and calcu- 
lating the dispersion of the best fit power law normalization 
and slope among the samples. 

Figure [2 and Table |2 show that the slope and evolution of 
the M^oQ-Tx relation are quite close to the self-similar model. 
There is a ~ 20% scatter in M around the mean relation and 
much of the scatter is due to unrelaxed clusters. Note also that 
the normalization of the M500 - Tx relation for unrelaxed sys- 
tems is somewhat biased with respect to that for the unrelaxed 
subsample. Unrelaxed clusters have somewhat lower temper- 
atures for a given mass. This may seem counter-intuitive at 
first, given that one can expect that shocks can boost the ICM 
temperature during mergers. However, in practice the effect of 
shocks is relatively small (e.g., O'Hara et al. 2006). The main 
source of the bias is that during advanced mergers the mass 
of the system already increased but only a fraction of the ki- 
netic energy of merging syste ms is converted into the th ermal 
energy of the ICM (see, e.g., Math iesen & Evrardi2001h . 

The Mg - M500 relation (Fig. |2j has a somewhat smaller 
scatter (a; 10 - 12%) around the best fit power law rela- 
tion than the M - Tx, but its slope is significantly different 
from the self-similar prediction — we find Msoo oc ^0.88^0.92 
compared to the expected M500 °^ Mg. This is due to the 
trend of gas fraction with cluster mass, /gas = Mg/M^oi) oc 
■^^500 °^ present for both the simulated clusters in our sam- 
ple (s ee iKrayteoy e^: al.,.2005) and for the observed clusters 
( Vikh linin et alJl2006h The normalization of the Mg - M500 
relation, on the other hand, evolves only weakly between 
z - 0.6 and z - Q (yet, statistically significant evolution 
is present; see Table IS, reflecting slow evolution of the gas 
fraction with time ( Kravtsov et al. 2005) . 

The Fx -M500 relation (Fig.|3jl has the smallest scatter, only 
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Fig. 1. — Relation between the X-ray spectral temperature, Tx, and total 
mass, Mjoo- is measured within the radial range [0.15 - IJrsoo. Separate 
symbols indicate relaxed and unrelaxed clusters, and also 2 = and z = 0.6 
samples. The dashed line shows the power law relation with the self-similar 
slope fit to the entire sample, and the dotted lines indicate 20% scatter. 



^ 5 - 7%. Note that this value of scatter includes clusters at 
both low and high-redshifts and both relaxed and unrelaxed 
systems. In fact the scatter in Yx - M500 for relaxed and un- 
relaxed systems is indistinguishable within the errors (see Ta- 
ble|3- Note also that the figures include points corresponding 
to the three projections of each cluster Figure|3]shows that the 
dispersion in the projected values of Yx for each given cluster 
is very small, which means that Yx is not very sensitive to the 
asphericity of clusters. Remarkably, the scatter of Yx - M500 
is as small as that in the three-dimensional Fsz - ^#500 relation 
(see Table|2j. This implies that Yx is in fact a better mass indi- 
cator than SZ flux, when additional possible sources of scatter 
for the real measurements of the SZ flux, such as projecting 
structures outside rsoo, are taken into account. 

The comparison of the mass proxies, clearly shows that Yx, 
the product of gas mass and X-ray spectral temperature, is by 
far the most robust and most self-similar mass indicator Why 
is the product better than its parts? Figure 0] shows that the 
answer lies in the anti-correlation of the residuals of temper- 
ature and gas mass from their respective relations with total 
cluster mass. We plot residuals from the best fit power laws 
with the slope value fixed to the self-similar value (i.e., the 
upper portion of Table|3 to illustrate both random scatter and 
systematic deviations from self-similar behavior. 

Figure |4] also shows that the clusters with temperatures 
lower than the mean temperature for a given total mass tend 
to have gas mass higher than the mean, and vice versa. Note 
also that there is some redshift evolution between z = and 
z - 0.6 — more clusters have negative deviations of temper- 
ature and positive deviations of measured gas mass at z = 0.6 
compared to z = 0. This redshift evolution is thus in the op- 
posite direction for the gas mass and temperature deviation. 
The measured Mc, systematically increases at higher z for a 
fixed total mass, because the clusters become less relaxed on 
average. For unrelaxed clusters, the ICM density distribution 
is non-uniform which results in over estimation of Mg from 
the X-ray data JMathiesen et al.ll999l) . Some of the decrease 
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Fig. 2. — Correlation between gas mass and total mass of the clusters. Both 
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the same as in Fig.0 The dotted lines indicate 15% scatter 
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Fig. 3. — Yx - Msoo correlation. The meaning of the symbols and lines is 
the same as in Fig.Q The dotted lines indicate 8% scatter 

of Mg at lower z may be due to continuing cooling of the ICM 
which decreases the mass of hot. X-ray emitting gas. 

The anti-correlation of residuals and opposite evolution 
with redshift for gas mass and temperature is the reason why 
the behavior of their product, on average, has smaller scatter 
and is closer to the self-similar expectation in both the slope 
and evolution. We discuss the origin of this behavior further 
in §|6lbelow. 

5. PRACTICAL ALGORITHM FOR ESTIMATING 
CLUSTER MASS USING Yx 

Suppose we have ix -M500 relation pre-calibrated by some 
external means: M500 - C E(zy Y^, and we would like to use 
it to estimate Msoq. The pre-calibration can be done using a 
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TABLE 2 

Power-law M500 = CX"^ best fit parameters and scatter 



sample" 


quantity'' 


'Wt.soo - Tx 


Mt,500 - Mg,5oo 


'Wt.soo - Yx 




all 2, 
all clusters 


logio C 

a 
scatter 


14.41 ±0.014 
1.500 
0.196 


14.35 ±0.010 
1.000 
0.144 


14.27 ± 0.006 
0.600 
0.076 


14.27 ± 0.006 
0.600 
0.070 


all 2, 
relaxed 


logio C 

a 
scatter 


14.36 ±0.016 
1.500 
0.137 


14.36 ±0.019 
1.000 
0.157 


14.26 ± 0.009 
0.600 
0.070 


14.26 ±0.010 
0.600 
0.080 


all 2, 
unrelaxed 


logio C 

a 
scatter 


14.44 ±0.018 
1.500 
0.191 


14.35 ±0.011 
1.000 
0.131 


14.28 ± 0.007 
0.600 
0.075 


14.27 ± 0.006 
0.600 
0.059 


2 = 0, 

all clusters 


logio C 

(X 
scatter 


14.39 ± 0.020 
1.500 
0.221 


14.37 ±0.013 
1.000 
0.131 


14.27 ± 0.007 
0.600 
0.068 


14.27 ± 0.008 
0.600 
0.078 



z = 0.6, logio C 14.43 ± 0.017 14.34 ± 0.015 14.27 ± 0.009 14.26 ± 0.007 
all clusters a 1.500 1.000 0.600 0.600 

scatter 0.164 0.151 0.084 0.059 



all 2, 
all clusters 


logio C 

(X 
scatter 


14.41 ± 0.014 
1.521 ±0.062 
0.195 


14.35 ± 0.008 
0.921 ± 0.023 
0.107 


14.27 ± 0.006 
0.581 ±0.009 
0.071 


14.27 ± 0.005 
0.585 ± 0.010 
0.067 


all 2, 
relaxed 


logio C 

(X 
scatter 


14.36 ±0.017 
1.533 ±0.103 
0.136 


14.36 ±0.015 
0.898 ±0.051 
0.115 


14.26 ± 0.008 
0.579 ±0.012 
0.053 


14.26 ± 0.009 
0.564 ±0.014 
0.058 


all 2, 
unrelaxed 


logio C 

a 
scatter 


14.44 ±0.018 
1.553 ±0.063 
0.186 


14.34 ±0.010 
0.931 ±0.029 
0.095 


14.28 ± 0.008 
0.589 ±0.010 
0.072 


14.27 ± 0.006 
0.600 ±0.010 
0.059 


2=0, 

all clusters 


logio C 

a 
scatter 


14.39 ±0.019 
1.524 ±0.070 
0.219 


14.37 ±0.012 
0.917 ±0.028 
0.090 


14.27 ± 0.007 
0.583 ±0.010 
0.064 


14.28 ± 0.008 
0.584 ±0.013 
0.075 


2 = 0.6, 
all clusters 


logio C 
scatter 


14.44 ±0.018 
1.590 ±0.086 
0.157 


14.31 ±0.009 
0.871 ±0.033 
0.077 


14.27 ± 0.008 
0.571 ±0.016 
0.075 


14.26 ± 0.007 
0.577 ± 0.012 
0.051 



Power law fits were performed for different subsets of our cluster sample, split in redshift 
and/or dynamical state. 



For each observable X (= Tx, Mg, Yx, Y^^), we fit power law relation of the form M500 = 

aX/Xo)"', with Xo = 3.0, 2 X lO" Mq, 4 x lO", and 5 X 10"^ Mpc^ for Tx, Mg, Yx, Fsz, 
respectively. 



well-observed sample of relaxed clusters or simulations. Note 
that the definition of Yx includes spectral temperature and gas 
mass within rsoo- In practice, however, when we want to use 
Yx to estimate Msoo we do not know the value of r^oo a pri- 
ori. This is not a problem for measurements of Tx, because 
the integrated spectral temperature is dominated by the signal 
from the inner, brighter regions of the clusters and is thus not 
sensitive to the exact choice of th e outer radius. For exa mple, 
we checked that for clusters in the lVikhlininetal] (l200g ) sam- 
ple, a fairly drastic change of the outer radius to O.Srsoo from 
the fiducial choice of r5oo adopted in our analysis, results in a 
systematic increase of the mean spectral temperature of only 
« 5%. The sensitivity to the exact choice of the inne r reg ion 
is also comparably small, as was discussed above in § 13.31 

The measurement of gas mass, however, will depend quite 
sensitively on the adopted outer radius, because the profile 
Mg{r) does not converge at large r. In practice, therefore, we 
have to solve for r^ix, simultaneously with estimating M500, 
which can be done with the following iterative algorithm. An 
initial value for Tx can be obtained by using the integrated 
temperature (including the center) to estimate an approximate 



value of rsoo through a crude mass-temperature relation, and 
then re-measuring the spectrum in the radial range (0.15 - 
l)r5oo. Then, using identity Msoo - 4/3 tt 500 pc(2) '"500' we 
can re-write the Yx - M relation as 

4/3 7r500p,(z)r5\o = C£(zf [Fx(r50o)]", (7) 

where Yx{r) = Tx Mg(r) and the constant C is from pre- 
calibrated Yx - M500 relation. Equation is used to solve 
for rsoo. A second iteration can then be done using Tx re- 
measured in the newly estimated radial range (0.15 - 1) r^oo, 
continuing the procedure until convergence. Once rjoo is de- 
termined, Msoo = 4/37r500pc(z)''5()o- 

It is useful to know how the observational calibration of 
the normalization constant C in equation0scales with the as- 
sumed value of the Hubble constant. Cluster masses deter- 
mined dynamically or through gravitational lensing scale as 
M oc This implies that the angular size that corresponds 
to rsoo is /i-independent. The gas mass derived from X-ray 
data for a fixed angular radius scales as Mg oc h^^/^. There- 
fore, the observed normalization of the Yx - M relation with 
the slope close to self-similar scales as C oc /z'/^. 
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Fig. 4. — Fractional deviations in temperature and gas mass for fixed A/toi 
relative to their respective best fit self-similar relations, M500 ^ and 
A/500 ^g- The fit includes all systems, both at z = (solid circles) and 
z = 0.6 (open circles). Note that the deviations for gas mass and tempera- 
ture are generally anti-correlated: clusters with large positive (negative) de- 
viations in Mg tend to have negative (positive) deviations in Tx- A similar 
anti-correlation exists in the trend with redshift (compare the distribution of 
points for z = and z = 0.6). 



6. DISCUSSION AND CONCLUSIONS 

In this paper we presented comparison of several X-ray 
proxies for the cluster mass — the spectral temperature and 
gas mass measured within rsoo, and the new proxy, Yx, de- 
fined as a simple product of Tx and Mg. Analogously to the 
integrated Sunyaev-Zel'dovich flux. Fx is related to the total 
thermal energy of the ICM. To test these mass proxies, we 
use a sample of clusters simulated in the concordance ACDM 
cosmology. The simulations achieve high spatial resolution 
and include radiative cooling, star formation, and other pro- 
cesses accompanying galaxy formation. We construct mock 
Chandra images for the simulated clusters and derive the X- 
ray proxies from a procedure essentially identical to that used 
in the real data analysis. 

The main result of this study is that Yx is a robust mass 
indicator with remarkably low scatter of only « 5 -7% in M500 
for fixed Yx, regardless of whether the clusters are relaxed or 
not. In addition, we show that the redshift evolution of the 
Yx - Msoo relation is close to the self-similar prediction given 
by equation |5l which makes this indicator a very attractive 
observable for studies of cluster mass function with the X-ray 
selected samples. 

The Tx-MjOQ relation has the largest scatter 20%), most 
of which is due to unrelaxed clusters. The unrelaxed clusters 
have temperatures biased low for a given mass. This is likely 
because during advanced, ongoing mergers, the mass of the 
system already increased but only a fraction of the kinetic en- 
ergy of merging systems is converted into the thermal energy 
of gas, due to incomplete relaxation. 

The Mg - M500 relation shows an intermediate level of scat- 
ter, «s 10-12%. This relation does not appear to be sensitive to 
mergers. It does, however, exhibit significant deviations from 
self-similarity in its slope, which is due to the dependence of 
gas fraction within r^oo on the cluster mass (iKravtsov et al.l 



l2005h. A simila r dependence exists for the observed clusters 
dVikhUnin et"an r2006) and we can thus expect similar trends 
in the Mg - M500 relation for real clusters. 

Generally, all the observable-mass relations we tested 
demonstrate a remarkable degree of regularity of galaxy clus- 
ters as a population. Tx, Mg, and Yx all exhibit correla- 
tions with M500 which are close to the expectation of the self- 
similar model, both in their slope and evolution with time, 
within the uncertainties provided by our sample. The only 
exception is the slope of the Mg - M500 relation. 

Given that our analysis relies on cosmological simulations, 
it is reasonable to ask whether the simulated clusters are re- 
alistic. Although simulations certainly do not reproduce all 
of the observed properties of clusters, especially in their core 
regions, the ICM properties outside the core in simulations 
and observations agree quite well. We illustrate this in Fig.|5l 
which shows that the M^ ^oo - Tx relations for simulated and 
observed clusters (the sample of relaxed clusters at z ^ Q from 
fVikhlinin et al. 2006). For this comparison we use only those 
clusters from simulations that appear regular and relaxed in 
their X-ray surface brightness image. Clearly, both simu- 
lated and observed cluste rs exhibit tig ht correlations between 
M^,50 and Tx (see also iMohret al.1 [19991 IVoevodkin eisll 
l2002^ which agree remarkably in their slope (Mg oc r'-^^) 
and normalization^. The normalizations derived from simu- 
lated and real clusters agree to ^ 10%, while slopes are indis- 
tiguishable and both deviate significantly from the expected 
self-similar value of L5. This is a consequence of significant 
trends in the gas fraction with cluster mass, M^oZ-^soo 
^500"°'^^ for bo th simulated (IKravtsov et all [20051) and ob- 
served clusters ( IVikhlinin et al.ll2006 l). The deviations from 
the self-similar model also manifest themselves in the absence 
of any noticeable evolution with redshift^. Interestingly, the 
real clusters show a similarly weak evolution in the Mg - T 
relation ( Vikhlinin et al. 2002). Figurel^shows that the likely 
explanation is that the clusters at z = 0.6 tend to be colder 
for the fixed Mtot but have higher estimated Mg 500 then their 
counterparts at z = because they are less relaxed. 

A similar level of agreement between the simulations and 
latest Chandra measurements exists also for the total mass 
vs. temperature relation, M500 - Tx- In fact, the normaliza- 
tion for our simula ted sample (Tab led agrees with the obser- 
vational results of 'Vikhli nin et alJ (|2006) to k 10%. This is 
a considerable improvement given that significant disagree- 
ments existed just several years ago (see § The resid- 
ual systematic 10% difference in the normalization is likely 
caused by non-thermal pressur e support from bulk gas mo- 
tions (Falte nbacher et alJ l2005t iRasia et alJ l2006t iLau et alJ 
2006), which is unaccounted for by the X-ray hydrostatic 
mass estimates. 

In Figurel6lwe compare the Yx - Mgoo relation for th e simu- 
lated c lusteri^ and for the Chandra sample of Vikhlini n et alJ 
('2OO6'). The observed clusters show a tight correlation with 
the slope close to the self-similar value. There is « 15% 
difference in normalization, likely explained also by neglect- 
ing the turbulent pressure support in the Chandra hydrostatic 
mass estimates. 

* For this comparison, the gas masses of simulated clusters are rescaled by 
a factor of 0.17/0.143 = 1.19 to reflect the dift'erence between the universal 
baryon f ractions adopted in the simulation and the value measured by the 
WMAP ( Snergel et al '200.3: .Tegmark et aU2 004i. 

^ Note that Mg in Fig|5]is not multiplied by the E(z) factor unlike the total 
mass in FigfTlandlsl 
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Fig. 5. — Relation betwen X-ray spectral temperature and gas mass for the 
relaxed subsample o f simulated clus t ers (ci rcles) and for a sample of relaxed 
Cliandra clusters of IVikhlinin et alj <200fl stars). Both gas mass and tem- 
perature are the quanities derived from analysis of real and mock X-ray data. 
The error bars in the Chandra measurements are comparable to the symbol 
size and are not shown for clarity. The gas masses for the simulated clusters 
are rescaled by a factor of 0.17/0. 143 = 1.19 to reflect the difl'erence between 
the universal baryon fractions adopted in the simulation and the value mea- 
sured by the WMAP. The dashed line shows the best fit power law relation 
with the slope 1.75. 

The excellent agreement of simulations and observations in 
terms of the relation between the two X-ray observables used 
to compute Yx (Mg - Tx) and a relatively good agreement in 
the Tx - MsoQ and Yx - M500 relations, gives us confidence 
that the results presented in this paper are sufficiently realis- 
tic. One can ask whether the real clusters show the same trend 
of decreasing scatter when Mg and Yx are used as mass indi- 
cators instead of Tx- Unfortunately, the existing data cannot 
answer this question because the mass measurement uncer- 
tainties for individual clusters are of order or larger than the 
expected scatter for relaxed clusters. 

Our results show that Yx is clearly most robust and most 
self-similar X-ray cluster mass indicator The biases exist- 
ing in mass estimates based on Mg and Tx anti-correlate both 
for a given redshift and in terms of evolutionary trends (see 
Figure |5j- This explains why their product, Yx, is a better 
mass indicator than Tx and Mg individually. The quality of 
Yx compares well to that for the actual three-dimension inte- 
gral of the ICM thermal energy (proportional to Y^z) in terms 
of its low scatter and self-similarity (see Table O. Yx may 
prove to be an even better mass proxy than Ysz, given that we 
use ideal 3D measurement of the latter while reproducing the 
actual data analysis for the former. Note also that Fsz is more 
sensitive to the outskirts of clusters, because it involves gas 
mass-weighted temperature (as opposed to the spectral tem- 
perature more sensitive to the inner regions), and thus should 
be more prone to projection effects. 

Note that Yx is also an attractive mass proxy from the data 
analysis point of view. First, it reduces statistical noise by 
combining the two independently measured quantities, Mg 
and Tx, into a single quantity. Consider for example how 
mass estimates are affected by Tx, the parameter that is the 
most difficult to measure and which is most affected by the dy- 
namical state of cluster. A 10% observational uncertainty or 
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Fig. 6. — Fx — !^sm relation for the z = sample of the simulated clus- 
ters (c ircles) and for a sample of relaxed Chandra clusters of Vikhlinin et afl 
KOOai stars). The gas masses for the simulated clusters are appropriately 
rescaled (see caption to Fig|5). The dot-dashed line shows the best fit power 
law relation for the simulated clusters with the slope fixed to the self-similar 
value of 3/5. The dashed line shows the same best fit power law, but with the 
normalization scaled down by 15%. 



a real deviation due to a merger in Tx translates into a ~ 15% 
mass uncertainty through the M - Tx relation and only 6% 
uncertainty through the Yx - M relation. Yx is also less sen- 
sitive to any errors in the absolute calibration of the X-ray 
telescope. For example, a mis-calibration of the low energy 
effective area typically translates into Tx and Mg measure- 
ment errors of the opposite sign, 6T/T ^ -26Mg/M„. These 
errors partially cancel mYx,dY/Y ^ 0.5 ^r/T, and are further 
reduced in the mass estimate: 5M/M = 0.6 6Y/Y * 0.3 dT/T. 
The error in mass is larger when other proxies are used, 
6M/M = 6Mg/Mg ^ 0.5 6T/T for the Mg - M relation and 
5M/M = 1.5 6T/f for the M - T relation. 

The robustness and low scatter make Yx an excellent mass 
indicator for observational measurements of cluster mass 
function at both z = and higher redshifts. The necessary 
data — an X-ray brightness profile and a wide-beam spectrum 
excluding the core — are easily obtained with sufficiently 
deep observations with Chandra, XMM-Newton, and Suzaku 
telescopes. The small scatter and simple, nearly self-similar 
evolution of the Yx - M relation hold promise for the self- 
calibration strategies for future large X-ray cluster surveys. 
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dation (NSF) under grants No. AST-0206216 and AST- 
0239759, by NASA through grant NAG5- 13274, and by the 
Kavli Institute for Cosmological Physics at the University 
of Chicago. AV is supported by the NASA grant NAG5- 
9217 and contract NAS8-39073. The cosmological simula- 
tions used in this study were performed on the IBM RS/6000 
SP4 system (copper) at the National Center for Supercom- 
puting Applications (NCSA). D.N. is supported by a Sher- 
man Fairchild postdoctoral fellowship at Caltech. We have 
made extensive use of the NASA Astrophysics Data System 
and arXiv . org preprint server. 



A NEW X-RAY CLUSTER MASS INDICATOR 



9 



REFERENCES 



Allen, S. W., Schmidt, R. W., Fabian, A. C, & Ebeling, H. 2003, MNRAS, 

342, 287 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Anninos, R & Norman, M. L. 1996, ApJ, 459, 12 
Amaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 
Borgani, S., Murante, G., Springel, V., Diaferio, A., Dolag, K., Moscardini, 

L., Tormen, G., Tornatore, L., & Tozzi, R 2004, MNRAS, 348, 1078 
Borgani, S., Rosati, P., Tozzi, P., Stanford, S. A., Eisenhardt, P. R., Lidman, 

C. , Holden, B., Delia Ceca, R., Norman, C., & Squires, G. 2001, ApJ, 561, 
13 

da Silva, A. C., Kay, S. T., Liddle, A. R., & Thomas, P A. 2004, MNRAS, 
348, 1401 

Dave, R., Katz, N.. & Weinberg, D. H. 2002, ApJ, 579, 23 

David, L. P, Slyz, A., Jones, C., Forman, W., Vrtilek, S. D., & Amaud, K. A. 

1993, ApJ, 412, 479 
Evrard, A. E., Metzler, C. A., & Navarro, J. F. 1996, ApJ, 469, 494 
Faltenbacher, A., Kravtsov, A. V., Nagai, D., & Gottlober, S. 2005, MNRAS, 

358, 139 

Ferland, G. J., Korista, K. T., Vemer, D. A., Ferguson, J. W., Kingdon, J. B., 

& Vemer, E. M. 1998, PASP, 110, 761 
Finoguenov, A., Reiprich, T. H., & Bohringer, H. 2001, A&A, 368, 749 
Gotflober, S., Klypin, A., & Kravtsov, A. V. 2001, ApJ, 546, 223 
Haardt, F & Madau, P 1996, ApJ, 461, 20 

Hallman, E. J., Mod, P M., Bums, J. O., & Norman, M. L. 2006, ApJ 

submitted (astro-ph/0509460) 
Henry, J. P. 2000, ApJ, 534, 565 
Henry, J. R & Amaud, K. A. 1991, ApJ, 372, 410 
Hu, W. 2003, Phys. Rev. D, 67, 081304 

Ikebe, Y., Reiprich, T. H., Bohringer, H., Tanaka, Y., & Kitayama, T. 2002, 

A&A, 383, 773 
Kaiser, N. 1986, MNRAS, 222, 323 
— . 1991, ApJ, 383, 104 
Kennicutt, R. C. 1998, ApJ, 498, 541 

Kravtsov, A. V. 1999, PhD thesis. New Mexico State University 

Kravtsov, A. V., Klypin, A., & Hoffman, Y. 2002, ApJ, 571, 563 

Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 

Lau, E., Kravtsov, A. V., & Nagai, D. 2006, ApJ in preparation 

Levine, E. S., Schulz, A. E., & White, M. 2002, ApJ, 577, 569 

Lewis, G. F, Babul, A., Katz, N., Quinn, T., Hemquist, L., & Weinberg, D. H. 

2000, ApJ, 536, 623 
Lima, M. & Hu, W. 2004, Phys. Rev. D, 70, 043504 
— . 2005, Phys. Rev. D, 72, 043006 
Majumdar, S. & Mohr, J. J. 2003, ApJ, 585, 603 
— . 2004, ApJ, 613, 41 
Markevitch, M. 1998, ApJ, 504, 27 

Markevitch, M., Bautz, M. W., Biller, B., Butt, Y, Edgar, R., Gaetz, T., 
Garmire, G., Grant, C. E., Green, P., Juda, M., Plucinsky, P. P., Schwartz, 

D. , Smith, R., Vikhlinin, A., Virani, S., WargeUn, B. J., & Wolk, S. 2003, 

ApJ, 583, 70 

Markevitch, M., Forman, W. R., Sarazin, C. L., & Vikhlinin, A. 1998, ApJ, 

503, 77 



Mathiesen, B., Evrard, A. E., & Mohr, J. J. 1999, ApJ, 520, L21 
Mathiesen, B. F. & Evrard, A. E. 2001, ApJ, 546, 100 
Mazzotta, P., Rasia, E., Moscardini, L., & Tormen, G. 2004, MNRAS, 354, 
10 

Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517, 627 
Motl, P M., Halhnan, E. J., Bums, J. O., & Norman, M. L. 2005, ApJ, 623, 
L63 

Muanwong, O., Thomas, P A., Kay, S. T., & Pearce, F. R. 2002, MNItAS, 

336, 527 

Nagai, D. 2006, ApJ submitted (astro-ph/()5 12208) 

Nevalainen, J., Markevitch, M., & Forman, W. 2000, ApJ, 532, 694 

O'Hara, T. B., Mohr, J. J., Bialek, J. J., & Evrard, A. E. 2006, ApJ in press 

(astro-ph/05 10064) 
Oukbir, J. & Blanchard, A. 1992, A&A, 262, L21 

Peebles, P. J. E. 1993, Principles of physical cosmology (Princeton Series in 

Physics, Princeton, NJ: Princeton University Press, — cl993) 
Pierpaoli, E., Borgani, S., Scott, D., & White, M. 2003, MNRAS, 342, 163 
Press, W. H., Teukolsky, S. A., Vettering, W. T., & Flannery, B. P 1992, 

Numerical Recipes in FORTRAN (Cambridge: Cambridge Univ. Press) 
Rasia, E., Ettori, S., Moscardini, L., Mazzotta, P., Borgani, S., Dolag, K., 

Tormen, G., Cheng, L. M., & Diaferio, A. 2006, ArXiv Astrophysics e- 

prints 

Rasia, E., Mazzotta, P., Borgani, S., Moscardini, L., Dolag, K., Tormen, G., 

Diaferio, A., & Murante, G. 2005, ApJ, 618, LI 
Reiprich, T. H. & Bohringer, H. 2002, ApJ, 567, 716 
Seljak, U. 2002, MNRAS, 337, 769 

Spergel, D. N., Verde, L., Peiris, H. V, Komatsu, E., Nolta, M. R., Bennett, 
C. L., Halpem, M., Hinshaw, G., Jarosik, N., Kogut, A., Limon, M., Meyer, 
S. S., Page, L., Tucker, G. S., Weiland, J. L., WoUack, E., & Wright, E. L. 
2003, ApJS, 148, 175 

Stanek, R., Evrard, A. E., Bohringer, H., Schuecker, P, & Nord, B. 2006, ApJ 
submitted (astro-ph/0602324) 

Tegmark, M., Strauss, M. A., Blanton, M. R., Abazajian, K., Dodelson, S., & 
et al. 2004, Phys. Rev. D, 69, 103501 

Vikhlinin, A. 2006, ApJ in press (astro-ph/0504098) 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, M., Miuray, 

S., & Van Speybroeck, L. 2006, ApJ in press (astro-ph/0507092) 
Vikhlinin, A., McNamara, B. R., Forman, W., Jones, C., Quintana, H., & 

Hornstrap, A. 1998, ApJ, 502, 558 
Vikhlinin. A.. VanSpeybroeck, L., Markevitch, M., Forman, W. R., & Grego, 

L. 2002, ApJ, 578, L107 
Vikhlinin, A., Voevodkin, A., MuUis, C. R., VanSpeybroeck, L., Quintana, H., 

McNamara, B. R., Gioia, I., Homstmp, A., Henry, J. P., Forman, W. R., & 

Jones, C. 2003, ApJ, 590, 15 
Voevodkin, A. & Vikhlinin, A. 2004, ApJ, 601, 610 

Voevodkin, A. A., Vikhlinin, A. A., & Pavlinsky, M. N. 2002, Astronomy 

Letters, 28, 366 
Voit, G. M. 2005, Reviews of Modem Physics, 77, 207 
Wang, S., Khoury, J., Haiman, Z., & May, M. 2004, Phys. Rev. D, 70, 123008 



